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Abstract. We show that the transition probabiUty of the Markoc chain 
{G{j, 1), . . . ,G(j,n))j>i, where the G{i,j)'s are certain directed last-passage 
times, is given by a determinant of a special form. An analogous formula has 
recently been obtained by Warren in a Brownian motion model. Furthermore 
we demonstrate that this formula leads to the Meixner ensemble when we 
compute the distribution function for G(m, n). We also obtain the Fredholm 
determinant representation of this distribution, where the kernel has a double 
contour integral representation. 



1. Introduction 

The starting point for the present paper are some nice resuhs from the interesting 
paper [18 by J. Warren. Let Bk{t), 1 < fc < n, be independent Brownian motions 
started at the origin and define Xk (t) , /c = 1 , . . . , n, recursively by 

(1.1) Xkit) = sup (Xk-iis) + Bk{t) - Bki-s)), 

0<s<t 

t>0. The muhi-dimensional Markov process X(i) = {Xi{t), . . . ,X„(t)) at a fixed 
time is closely related to the largest eigenvalues of succesive principal submatrices 
of a GUE matrix. In fact, let H — (/iij)i<ij<n, be an n x n GUE matrix, i.e. 
distributed according to the probability measure exp{—Tr H^)dH on the space 
of n X n Hermitian matrices, and let Hk = {hij)i<ij<k, I < k < n, he the principal 
submatrices. Then, if XmnKiM) denotes the largest eigenvalue of the Hermitian 
matrix M, we have X{l/2) = ( Amax ),..., Ainax(i?n)) in distribution, [1], [6], 
[18j . Furthermore, there is a nice formula for the transition function of the Markov 
process X(t), [18], 

(1.2) P[X(i) = y I X(s) = x]= det(Z?J-Vt_,(2/, - a;,))i<,,j<„, 

if a;i < • • • < Xn, yi < • • • ^ Vn- Here D denotes ordinary differentiation, D^^ is 
anti-derivation, 

(1.3) 

and (jitix) = (27r<)^^/^ exp(— a;^/2t) is the transition density for Brownian motion. 

Let JGUE(n)('7) be the distribution function for the largest eigenvalue of an n x ti 
Hermitian matrix H from GUE. It follows easily from (|1.2p that 



(1.4) FG^Ein)iv) = det(i?^-^+Vi/2(^))i<»,,<«. 
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This formula, given in [T^ , can also be obtained directly from the GUE eigenvalue 
measure, see proposition 12.31 below. 

We will show in this paper, starting from definitions, that we have analogous 
formulas for the vector G{i) = {G{i, 1), . . . , G{i, n)), i > of certain last-passage 
times defined as follows, [?]■ Let w{i,j), G Z'^, be independent geometric 

random variables with parameter g, < g < 1, 

(1.5) V[wit,j) = k] = {l-q)q'' 
and define 

(1.6) G{m,n) ~ max \ w{i,j), 

where the maximum is over all up/right paths from (1, 1) to (m, n). It is clear that 
these random variables satisfy the recursion relation 

(1.7) G{m, n) = max(G'(TO — 1, n), G{'m, n ~ 1)) + wlm, n), 

where G(0, n) — G{n, 0) = for n> \. If we use this recursion relation repeatedly 
we see that 

G{m,n) = max (G'(j, n — 1) + > w{i,n) — y w{i,n)), 

l<j<n ^ — ^ ^ — ^ 
i=l i=l 

which looks like a discrete version of (jl.ip . From this it is reasonable to expect 
that there should be a formula for the transition function for the Markov chain 
(G(i))i>o similar to (jl.2p . This is indeed the case in a very natural way where 
the differentiation operator is replaced by a finite difference operator, see theorem 
12.11 This will imply a formula similar to (jl.3p for the distribution function of 
G(m, n), see theorem 12.21 We will also show how we can go from this formula 
(12. 3p to the known expressions [13], [S], for this distribution function in terms of 
the Meixner ensemble, f7], and as a Fredholm determinant with a double contour 
integral expression for the kernel, [14] , [9] . The Fredholm determinant formula has 
the advantage that it is much better suited for computation of asymptotics. The 
argument in this paper leading to (|2.8p gives an alternative approach, starting from 
the definitions, to this formula. 

Results related to the transition probability (|2.2p for (G(i))i>o, but in the case of 
w{i,j) exponentially distributed, go back to the work of G. Schiitz, 16:, where the 
totally asymmetric exclusion process (TASEP) is studied using the Behte ansatz, 
see e.g. [7] for a discussion of the relation to G(m, n). This is not exactly the same 
Markov chain, but the results of Schiitz can also be used to derive the expression for 
the distribution function for G{m, n) in terms of the Laguerre ensemble, see |15| and 
also [13] . In [T5] the case of geometric random variables is also considered and the 
formula for the distribution function for G(m, n) in terms of the Meixner ensemble 
derived. This is based on results from [5] on a discrete TASEP-type model. More 
general formulas for the asymmetric exclusion process (ASEP) have been proved 
recently in [T^. Results generalizing the formula (dHj) to discrete models has also 
been given independently in [5 . 

2. Results 

For z G Z we define w{x) = (1 — q)q^H{x), < q < 1, where H is the Heaviside 
function, H{x) = if a; < and H{x) = 1 if x > 0. The m-fold convolution of w 
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with itself is then the negative binomial distribution, 

(2.1) w^{x) = {l-qr(''^"''^)q^H{x), 



as is not difficult to see using generating functions. 

For a function / : Z ^ C we denote by A the usual finite difference operator, 
A/(a;) = fix + 1) - fix). We also set 

{A-'f)ix)= ^ f{y) 

y=-oo 

provided the sum is convergent, and A°/ ~ f. Note that A(A^^/) = A^^(A/) = 
/• 

Let Wn ^ {x E ; xi < ■ ■ ■ < x„}, and let G(i), i > 0, be the Markov chain 
defined in the introduction. 

Theorem 2.1. If x,y g Wn, then for m> £, 

(2.2) P[G(m) = y \ G{£) = x] ^ det(A^-'u;„,_f (y, - a;,))i<»j<n. 

We postpone the proof to section [3] The proof first shows (|2.2p in the case 
m = £ + 1 using (|1.7p and an induction argument in n, and then establishes a 
convolution formula for the determinants involved using the generalized Cauchy- 
Binet identity ([23]) . 

Taking G(0) = it is not difficult to show, see section [3l that we have the 
following consequence of (|2.2p . which is analogous to (|1.4p . 

Theorem 2.2. For any rj E N, m> n> 1, 

(2.3) P[G(m, 7i) < 77] - dct{A^-'~^Wmiv + l))i<^j<n- 

As stated in the introduction it is possible to relate the expression in the right 
hand side of (|1.3p directly to the expression for the distribution function coming 
from the GUE eigenvalue measure. Let 

A„(x) = det(a;^"^)i<ij<„ 

be the Vandermonde determinant. 

Proposition 2.3. We have the following identity for any 77 € M, n > 1, 

(2.4) det(Z?^-^+Vi/2(77))i<^,,<„ = ^ / Anixf J] e~"?rf"x, 

J(-oo,J7]" ^^-^ 

where Zn is the appropriate normalization constant. 

Proof. Let Hj, j > 0, be the standard Hermite polynomials. Then Hjix) = 2^pjix), 
where pj is a monic polynomial, and we have Rodrigues' formula 

D^e-^' = {-iyHj{x)e-''\ 

Hence, 

D'-'-'^i/2iv) = D-\D=-'(l),/,)iv) ^ f ^^^^^^D^-'<P^,2ix)dx 
= / {x - rir-'p,ix)e-^ dx. 
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Using row and column operations we obtain 

det(D^-'-Vi/2(f]))i<»j<n = 11 -^det( / x'-^x^-^e-^ dx))l<^J<n 

j=0 ■ V J —oo 

1 r " 

J(-oo,J7]" ^-^j 

with Z„ = 2^"("^i)/27r"/^ rijCo^ i'- ^^"^ ^^^^ equality we have used the generalized 
Cauchy-Binet identity, 

(2.5) det{( (t>,{x)%l]j{x)d^{x)) = —^ j det(0j(xj)) det(V'i(a;j)) Jl d^(a;j) 



where all determinants are n x n. □ 

We have a similar identity relating the right hand side of (|2.3p to the Meixner 
ensemble. The proof is a little more involved and we postpone it to section [31 

Proposition 2.4. For any rj ^ N, m > n > 1, 
(2.6) 

detiA^-^-'w^r, + l))i<,,<„ = Yl ^"(^)' n f"^' ""V^^- 

' 0<Ki<r(+n— 1 j = l ^ J ' 

For asymptotic analysis it is more useful to have a representation of the distribu- 
tion function for G(to, n) as a Fredholm determinant with an appropriate kernel. It 
is possible to go to such a formula using Meixner polynomials and a standard ran- 
dom matrix theory computation as was done in f7]. There is also another formula 
for the kernel as a double contour integral which can be obtianed from the Schur 
measure, see [13], [S] or [TU]. It is actually possible to go directly to a Fredholm 
determinant formula with a double contour integral formula for the kernel starting 
from the expression in the right hand side of (|2.3p . 

Let 7r denote a circle centered at the origin with radius r > 0. Let 1 < r2 < 
ri < l/q and define 

, , , , 1 f dz f dw w (1 - gz)™(l - wY 



(27ri)2 z J^^^ w w — z wy+^^ (1 — z)"(l — gw)" 

for X, y G Z. 

Proposition 2.5. For any integer rj > 0, m > n > 1, 

(2.8) V[G{m,n) < ?]] ^ det{I - K,n,n)e^{,i+i,rj+2, ...})■ 

The proof will be given in the next section. 

Remark 2.6. It would be interesting to understand the joint distribution of G{mi, rii), 
i = 1, . . . ,p in order to understand the fluctuations of the "last-passage times sur- 
face", 3 {m,n) {G{m,n). If (mi,ni), . . . , (mp^rip) form a right/down path 
then the joint distribution of G(mi, ni), . . . , G{mp, rip) can be expressed as a Fred- 
holm determinant, see [TU], [3], and it is possible to investigate the asymptotic 
fluctuations. However, the asymptotic correlation between for example G(m, m) 
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and G{n, n), m < n, is not known and their is no nice expression for their joint dis- 
tribution. Using ()2.2|) we can write down an expression for their joint distribution, 
which was one of the motivations for the present work. We have 

P[G(m,m) < 77i,G(n,n) < 772] 
(2.9) 

= ^ ^ det(A-'~'i«„(xj))i<ij<„ det(A-'~'w„_„(2/j - a;i))i<ij<r 

X^^n, a;m<?7i y^Wn. ynl£V2 

However, we have not been able to rewrite this in a form useful for asymptotic 
computations. 

Remark 2.7. The case when the j)'s are exponential random variables can 
be treated in a completely analogous way or by taking the appropriate limit of the 
formula above, q ~ 1 — a/L, G{m, n) G{m, n)/L and I 00. 

Remark 2.8. Random permutations can be obtained as a limit of the above model, 
q = a/n^, n —> 00, 8 . The random variable G{n,n) then converges to L{a) the 
Poissonized version of the length ij^ of a longest increasing subsequence of a random 
permutation from Sn- We can take this limit in the formulas (|2.7[) and (|2.8[) and 
this leads to a formula for P[L(a) < 77] as a Fredholm determinant involving the 
discrete Bessel kernel, [8], [2]. Hence, we obtain a new proof of this result which 
does not inolve some form of the RSK-correspondence. 



3. Proofs 

3.1. Proofs of theorems 12.11 and 12.21 To prove theorem 2.1 we first consider 
the case m — £ = 1. The transition function from G(^) to G(m) is, by (|1.7p . 

TL 

(3.1) P[G{£ + 1) - y I G(£) = = [] w{yk - max(xfc, yfc_i)), 

where we have set yo = and x,y G Wn- Note also that it is clear from (|1.7[) that 
(G(i))i>o is a Markov chain. The right hand side can be written as a determinant 
by the following lemma. 

Lemma 3.1. Ifx,yG Wn, then 

n 

(3.2) Y]_ '^iVk - niax(a:fc,yfe_i)) = det(A^~*ii;(2/j - Xi))i<ij<n. 

Proof. We use induction with respect to n. the claim is trivial for n = 1. Assume 
that it is true up to n — 1. Expand the determinant in p.2p along the last row, 

dct(A^"*w(?/j - Xi))i<ij<„ = 

n-2 

^(-l)^-+"A'^-"u;(2/fc - Xn)det{A^-'w{y, - x,))i<,,,<„_i,,^fc 
fc=i 

- A^^w(j/„_i - Xn) det(det(A-'"'w(j/j - a;i))i<ij<„_ij^„_i 

(3.3) +uj(j/„ - 2;„)det(A^"*u;(?/j - Xi))i<ij<n-i. 
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We will first show that each term in the sum from k ~ 1 to 71 — 2 in the right 
hand side of (13. 3p is zero. Let Aj, denote the difference operator with respect to 
the variable y. The fact that A(A~^i«) = w, then gives 

det{A^~''w{yj - a;j))i<j,i<n-i,i/*; ^Vk+i ■ ■ ■ ^y^ det{A^'''w{yj - Xi))i<ij<n~i, 

where yj = if 1 < J < fc — 1, yj = yj+i if fc < n — 1. By the induction asumption 
this equals 

k-l n-l 

Aj,^^^ . . . Aj^„ Y[w{yj-max{xj,yj_i))w{yk+i-raax{xk,yk-i)) J| w{yj+i-ma.x{xj,yj)). 

j = l j=k+l 

If yk < Xn, then A''~"w{yk — Xn) = since A~^'w{x) = if a: < j- In this 
case the fc 'th term in the sum in (|3.3p is zero. Assume that y^ > Xn- Then 
?;„>•••> yfc > and we obtain 

det{A^'^w{yj - Xi))i<tj<n-i,j^k 

k-l ri-1 

= n^fe-max(:r„y,_i))(l-g)g— (--^-i)A,,^,...A,„ [] (1 - . 

j=l J=k+1 

But, 

n-l 

Ay,^, . . . A,„ n (1 - ih"'^'-'' = (1 - qr-'^'Ay,^, . . . Ay„qy- = 
j = k + l 

since fc + 1 < n. Hence, each term in the sum from fc = 1 to n — 2 in (|3.3p is zero 
and we obtain 

det(A-'~*w(?/j - Xi))i<tj<n 

= -A"^w(?;„_i - a:„) det(A-'"*w(?;j - Xi))i<ij<„_i j^„_i 
(3.4) + w{yn - Xn) det(A^"*w(yj - a;i))i<i j<„_i . 

If J/n-i < 2;„, then A~^i(;(y„_i — a;„) = and the right hand side of p.4p is 

n-l 

w(y„ ~ Xn) Y[ w{yj - max{xj,yj^i)), 
i=i 

by the induction assumption. Since yn-i < Xn, w{yn — Xn) = — max(a;„, y-n-i) 
and we get exactly the right hand side of (|3.2p . 

Assume now that > Xn- Note that A^w{x) = (q — l)A^^-'^w(a;) if j > 1 and 
a; > 0. The determinant 

det(A^"'w(yj - a:j))i<i<„_i j^„_i 

has A'^~'^w{yn — xi) = {q — l)A"^^^'w(y„ — Xi), 1 < z < n, in the last column, 
since ?/„ > J/n-i > > ■ ■ • > a^i- By the induction assumption this determinant 
equals 

ji-2 

(g- 1) J]^ w{yj - max{xj,yj^i))w{yn - max(a:„_i, ?;„_2))- 
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Thus, the right hand side of (|3.4p equals 

n-2 

- (1 - q){q'^"-^~''" ~ 1) n ^'^y^ ~ max(a;j,yj_i))w(y„ - max(a;„_i, y„^2)) 
(3.5) 

i=i 

where we have used the fact that (A~-^w)(y„_i — a;„)(g — 1) = (1 — — 1) 

for i/n-i — a^n > 0. since y„ > > a^n 2i ■ ■ ■ 2i the expression in p.5p can be 
written 

n-2 

X n ^'^y^ ^ max(x-,-,?/-,_i)) 

n — 1 n 

= (1 - q)qy"~y"-^ Y[ "^(yj ~ max(a;j,?/j_i)) = J| - max{xj,yj-i)). 

□ 

Theorem 12 . 1 1 foUows from lemma [XT] and the following convolution type formula 
for determinants of the form we have. 

Lemma 3.2. Assume that /, g : Z — )■ C are such that f{x) — g{x) — if x < M 
for some M G Z. Then, 

(3.6) det(A^-7(2;, -a;,))det(A^-'5(z, -y,)) = <lei{{A^-' f * g)(z, - x,)), 

where all determinants are n x n. 

Proof. We will make use of the following summation by parts formula 
(3.7) 

b b 

Au{y~x)v{z—y) = u{y—x)Av{z—y)+u{b+l—x)v{z—b)—u{a—x)v{z+l—a). 

The first step is to show that 

J2 det(AJ"-7(y, - x,)) det{A^-' g{z, - y,)) 
y&W„ 

(3.8) = J2 det(Ai-7(2;, - x^,)) det(A'-ig(z, - y,)) 

by repeated summation by parts. The left hand side of (j3.8p can be written 

^ A,„det(Ai-7(yi-:r,)...A"-i-7(2/n-i-a:,) A"-'-' f{y„ - x,)) 
y£W„ 

X det(A^-i6r(z, - yi) . . . A^-"g(z, - y„))- 

Here we have used the summation by parts formula (|3.7p to sum y„ between y,i_i 
and oo. The terms coming from u{b + 1 — x)v{z — fe) — u(a — x)v{z + 1 — a) in 
(|3.7p with b ^ oo and a = y„_i give since A'~"(ji(zi — 6) = if 6 is large enough 



8 K. JOHANSSON 

and the other term gives rise to a determinant with two equal columns containing 
A"-i-V(2/„-i - x^). The result is 

J2 det(Ai-7(2/i-x,)...A"-i-7(2/»-i-a;.) A^'-'-' fiy,, - x,)) 
y&W„ 

X det(A'-i5(z, - yi) . . . A'-"+\g(z. - Un-i) A'^'^+'giz, - y„)). 
We can now repeat this procedure with j/n-i, yn-2, ■ ■ ■ ,y2, which gives 
J2 det(Ai-7(yi - A'-'f{y2-x,) A^'^ /{y^ - x,) . . . A"-i-7(y„ - a;.)) 

xdet(A*-ig(^,-yi) A'-'g{z,-y2) A'~^g{z, - y^) . . . A^-"+ig(z, - y„)). 

Again we repeat the summation by parts procedure with j/„,...,y3, then with 
yn, ■ ■ ■ ,y4 and so on until we get the right hand side of (|3.8p . 

Next, we apply the generalized Cauchy-Binet identity p.Sp to the right hand 
side of dMl). This gives 

det(^ Ai-7(2/ " x.)A^-'giz, - y)). 

To prove the lemma it remains to show that 

^ Ai-7(y - x)A^-^g{z ~ y) = A^-\f * g)iz - x). 

If we set h{x) — H{x — 1), then A~^/(a;) = /i * f{x) and hence 

^'-'fiy ~ x)A^~'g{z -y) = A^\h<^-'^ * /) * giz - x) 

= ^-\A'-\f * g)){z -x) = A^-\f *g){z~x). 
Here h*^ denotes the j-fold convolution of h with itself. □ 
To prove theorem 12.21 we note taht by theorem 12. II 

P[G(m,n) < 77] = det{A^-'w^{xj)) 

= E E det(A^-^«;™(a;7). 

Now, 

J2 det{Ai~'wm{xj)) 

Xn—Xn-i 

= det{A^-'Wrn{xi) . . . A^-^-'w^{Xn-l) A^'^-^Wmiv + 1) - A^''^-' W^{Xn-l)) 

= det{A^-'wrn{xi) . . . A"-i-^i(;„(a;„„i) A"-i~V,„(77 + 1)). 
Repeated use of this argument proves theorem 12.21 
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3.2. Proofs of propositions [2741 and [2T5l The generating function for w{x) is 

^ ^ ' l-qz 
and lience, since i/7,„ is the m-fold convolution of w with itself, 



where the radius r of the circle 7^, centered at the origin, satisfies < r < l/q. 
It follows that 

(3.9) /^^w„,{x) = / ^l—lt^dz, 

^ ' ^ ' 2TTi (1 - 9z)™z^+'=+i ' 

for fc > (and also for fc < if r < 1). 

As noted above, if h{x) — H{x — 1), then A^^/(a;) — h * f{x) and hence 
ISr^ J{x) — [h*^ * f){^)- Using e.g. generating functions it is not difficult to 
see that 

(3.10) h*\x)^^-^^^^H{x-k), 

where y^^^ — y{y — 1) . . . (y — k + 1) is the factorial power. We can write 

Note that - y) = if x - y < i by p.lOp and hence h"{x-y) = if y > x - 1 
for any z > 1. We obtain 

(3.11) A^— iu;„(r; + l)^ ^ (!^_ll!!_!lA^-i«;„(y). 

y— — oo ^ ^' 

Fix L > n — 1. By (j3.1ip and some row operations we find 

^ {y_ + Ly-' 



detiA^-'-'wmiV + l))l<^,l<n = det( J2 y„ L hiy-'A^~'wrniy))l<^,j<n 

y — ~-oo ^ ' 

(3.12) = det( Iy±i^)(-i).'--iA^-iz«™(y))i<,,,<„, 

since it follows from (|3.9p that A^~^Wm^y) = if y < — n for \ < 3 <n. 

If we choose L = n — 1 it follows from theorem 12.21 (|3.12p and the generalized 
Cauchy-Binet identity (12. 5p that 
(3.13) 

n— 1 ^ jj+ri — 1 

P[G(TO,n) < ^] = n -j H A„(y)det((-l)^-iA^~iu;™(y,-n+l))i<,,,<„. 
If < fc < n, m > n, then 

m— 1 

(3.14) (-l)'=A'=u;,„(y) = Sfe(y) J] (y + ^)g'^H(2/ + n - 1), 

where Sfe is a polynomial of degree fc. To see this we can use p.9p . We see that 
the integral in (|3.9p is zero ifa;<— (fc + 1), and in particular if 2; < — n for all fc, 
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< k < n. If we make the change of variables z — > l/z and assume that x > 1 — m, 
we find 

This is a polynomial of degree m — 1 in y times , and this polynomial has zeros 
at —{k + 1), ... ,1 — m. Consequently p.l4p follows. 

Furthermore we have the following determinantal identity. Let pj, j = 0, . . . , n—1 
be polynomials of degree j and A2, . . . , An-i constants. Then there is a constant 
B such that 

n 

(3.15) det{pj_i{xi) Y\_ {xi + Ak))i<ij<n = BAn{x). 

k=j+i 

This is not hard to see. Choose Crj so that 

n n 
^X^'^'^Crj = Pj-l{x) Y[ (X + Ak). 
r—1 k—j+1 

Then, the left hand side of p.lSp is 

n 

detC^xl~^Crj)i<ij<n = An{x)detC, 
and we have proved (|3.15p with B = det C. Actually, according to [T2] we have 

n 

B = l[{-iy-'pj-i{-Aj), 
j=i 

but we will not need this result. 

By (IXT^ . (jXTil) and ((gl^ we obtain 

rj+n—l n m — 1 

P[G(m, n) < 77] = C,n,n E ^"(2^)' n n +k + l-n) 

yi,...,yn—0 j — 1 k—n 

77+71— 1 



' vi,---,yn=o j=i ^ 



Uj + m — n 



for we let 77 —> 00 we see that Zj^ ^ must be exactly 

the normalization constant in the Meixner ensemble. This proves proposition 12.41 
We now turn to the proof of proposition 12.51 Write K = m — n + 1 and define, 

for < j < n, X E Z, 

ai(x) — / z — — - — az 

' A„ (z- 1)^+1 



and 
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where 1 < r2 < ri < 1/q. If a; > 0, < j < n, then 

(q- iy+^ 

(3.16) — Pj{x ^ 

where pj is a monic polynomial of degree 1. This follows from the computation 
9-1 /■ .+„-i(g2-lP+^-^ 



9-1 



27ri 



"'7^2 r=0 ^ ^ ^ 



We see that this is a polynomial of degree j in a; with leading coefficient 
It follows from theorem 12.21 and p.l2p with L = —n that 

CV+n \ 



where, by p.9p . 

(3.18) i-iyA^w^iy-n)^^^—^ f J^—l)L^_^dz. 



Set 



Then, 

(3.19) E (gTi/^^+z^ ^'^^^^ " = ^'^y^- 

To show this it is sufficient to show that 

-^ (i-g)"'-^-^(-i)^-+^(z-i)^- [z-iy 

by (|3.18p and the definition of b^. The identity (|3.20p can be rewritten as 

n-l 

(3.21) ^(1 - -)^(1 - q)"-^-^-^(-iy+'<c,, = (-1)^(1 - qzr-'-\l - 
Set w = 1 — 1/z. Then (|3.2ip is equivalent to 

n-l 

^ w^l - (7)"-J-i(-l)^Cjf = (1 - g - u')"-^-iw^(-l)^. 
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By the binomial theorem 



n-1 



n— 1— r 



and we have proved ([XTQ]) . Note that det(c-,^) = 1 and hence by ([3?T6l) . ([3l7l) and 
P[G(m, n) < 77] = det ^ ^' ( " 

/v+n n-1 ^_^y 



det ( E E (^q_iy+K ^^^rniy - n)cji 



yy=o 



0<i,e<n 



(3.22) =det E«^(y)^j(j') 

We will now use the fact that 

00 

(3.23) ^aj{y)bk{y) = 6jk- 

v=o 

To prove this we use the definitions of aj and bk, 

X=0 ^ ' ''■2 'i-l \£C = 



(qz - iy+^-\w - 1)'= 
(z-l)J+i(gu;-l)'=+^ 



(7-1 /■ , f , w 1 (gz- l)^+^'-Vw- 1) 
' az aw 



K ■ 



(27ri)2 7^^^ J^^^ zw-z{z-iy+^{qw~l)''+ 

Since 1 < r2 < ri < 1/q we see that w = z is the only pole in the w-integral and 
hence by the residue theorem this equals 

27ri 



{qz-iy-''-Hz-l)''~^~^dz. 

J'1r2 

li j < k the integral is zero by Cauchy's theorem, li j > k we make the change of 
variables z 1/z and we see again that the integral is zero. When j — k we get 



g-1 
27ri 

Using (IX^ in ((X^ we get 



/ {qz-l)-^{z-l)-^dz^l. 



^[G{m,n) < 77] = det(% - E ai(y)^j(y))o<ij<n 

= det(/ - Km,n)p({-n+l.q+2.... })7 
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where 

q-l 



ji-l 

E 

J=0 



1 

27ri)2 
1 

27ri)2 



aj(a; + n)bj{y + n) 



— / dw 



7ri 



i^y+n (^qui- l)^(z- 1) ^ L(-2 - 1)('7«^ - 1) 



— / dw 



dz 
z 

dz 

z 



K-1 



K-1 



(1 - w)"(l - qzY' 



dw (1 - qz)^-^ (1 - - qz) 



(1 - z)"(l - qwy 
qz) 



w - 



w 



w wy+"- (1 - qw)^-^ (1 - z)"(l - qw)' 



w - 



Here we have used the fact that by Cauchy's theorem the — 1 term in the third hne 
gives a zero contribution in the z-integral since ri > r2- This proves proposition 



Remark 3.3. The Meixner ensemble is related to Meixner polynomials since these 
polynomials are orthogonal on N with respect to the negative binomial weight. 
These polynomials are not used explicitely in the computations above but figure in 
the background. This can be seen from the Rodrigues' formula, 



P 



im,q) 



{x)q- 



3+x 



m — 1 



TO — 1 
X 



llix-k) 



q I \ [x 

k=0 



and the integral representation 



(x) 



J 



(1 — z/q)^ dz 



' ' 27ri J^^ (1 - z)^+^ zi+i 



with < r < 1, where Pj''^ {x) are the standard Meixner polynomials, [TT] . 

Acknowledgement: I thank Jon Warren for a discussion and for sending me a 
preliminary version of [S] . 
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